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Magnetic quantum tunneling in a single molecule magnet (SMM) has been studied in experiments 
that combine microwave spectroscopy with high sensitivity magnetic measurements. By monitoring 
spin-state populations in the presence of microwave magnetic fields, the energy splittings between 
low lying high spin superposition states of SMM N14 (S = 4) have been measured. Absorption 
linewidths give an upper bound on the rate of decoherence. Pulsed microwave experiments provide 
a direct measure of the spin-lattice relaxation time, which is found to be remarkably long (~sec) 
and to increase with the energy splitting. 

PACS numbers: 75.45.+j, 75.50.Tt, 75.60.Lr 



Magnetic quantum tunneling (MQT) in single 
molecule magnets enables the creation of coherent high 
spin superposition states, which are both of great fun- 
damental interest and essential for the use of SMMs in 
quantum computing 0|. The main focus of experi- 
ments conducted to date, however, have been on incoher- 
ent MQT, where these superposition states are subject 
to rapid decay into their classical counterparts: spin-up 
and spin-down states 0.Q- Decoherence generally occurs 
when discrete levels are coupled to an environment with 
many degrees of freedom, such as the modes of a lattice 
(phononsjk the electromagnetic field (photons) or nuclear 
spins 0, IE • Here we investigate coherent MQT in a 
SMM in which the tunneling rate is faster than the rate 
of decoherence and at a temperature at which tunneling 
occurs only between the lowest lying spin states. Our ex- 
periments combine microwave excitations with high sen- 
sitivity magnetic measurements and enable continuous 
monitoring of spin state populations during the applica- 
tion of microwave radiation. Pulsed microwave exper- 
iments provide a direct measure of the spin-lattice re- 
laxation time. Further, the linewidths associated with 
microwave absorption provide an upper bound on the 
decoherence rate of a high spin superposition state in a 
SMM. 

We have chosen to study [Ni(hmp)(t-BuEtOH)Cl] 4 , 
henceforth referred to as N14, because this is a particu- 
larly clean SMM, with no nuclear spins on the transition 
metal sites 8]. The molecule core consists of four Ni 7/ 
(Spin 1) magnetic ions and oxygen atoms at alternating 
corners of a distorted cube, with S4 site symmetry. Fer- 
romagnetic exchange interactions between the Ni 11 ions 
lead to an S = 4 ground state at low temperature, as de- 
termined by magnetic susceptibility measurements and 
high frequency Electron Paramagnetic Resonance (EPR) 
8j . Magnetic hysteresis is observed below a characteris- 
tic blocking temperature (1 K) and is associated with the 
presence of a uniaxial magnetic anisotropy that favors 
high spin projections on the easy axis of the molecule, 
which we denote the z-axis. This uniaxial anisotropy 
leads to large energy barrier (DS 2 ~ 12 K) to magneti- 



zation reversal (see Fig. 1A). 

MQT is characterized by the spin Hamiltonian: 

H = -DS 2 Z - lib~S -g-H (1) 

where the first term is the uniaxial anisotropy and the 
second term is the Zeeman energy associated with the 
interaction between the spin and magnetic field. An ex- 
ternal magnetic field applied along the easy axis of the 
molecules H z , tilts the double potential well favoring 
those spin projections aligned with the field. Note that in 
zero magnetic field the \up) and \down) spin-projections 
have nearly the same energy and MQT is possible. Im- 
portantly, in this case, a magnetic field transverse to the 
anisotropy axis (in the x — y plane) lifts the degeneracy 
of these states by an energy A, the tunnel splitting, and 
leads to states that are coherent superpositions of the 
original \up) and \down) spin-projections. The frequency 
of MQT between \up) and \down) states is proportional 
to the tunnel splitting (/ = A/h). Fig. 1C shows the 
energy of the lowest lying states as a function of longitu- 
dinal magnetic field. At zero field the states are super- 
positions of opposite spin-projections 

\ S ) = Tad^) + \ down )) (0 s 
\A) = ^(\up)-\down)) [Z) 

where 

M = 7f ^m( fl ™ + t> m )\m) 

\down) = -^J2 m ( a m - b m )\m) 

The value of the coefficients a m and b m depend on the 
applied field. Only the coefficients with \m\ = 4 are non- 
zero at zero field and the two lowest levels are \S) — 
i(| + 4) + I - 4)) and \A) = (| + 4) - | - 4)). In 

the presence of a transverse magnetic field the \up) and 
\down) states are tilted away from the z-axis but are 
still separated by a large angle (see Figs. IB and 1C), 
when this field is less than the anisotropy field H a = 
2DS/(g^ B ) = 4.5 T. 
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FIG. 1: (Color on-line) A) Magnetic energy levels of Nu, 
with different z-axis spin projections labeled. B) The \up) 
and \down) state probability amplitudes versus z-axis projec- 
tion for different transverse magnetic fields. C) Energy of the 
lowest lying levels in the vicinity of zero field. D) Magneti- 
zation of lowest (black) and first excited states (gray) as a 
function of the longitudinal field for two different transverse 
fields. E) Experimental configuration. 

The interaction of the spin with the environment lim- 
its the coherence time of the superposition. A transverse 
field can be used to increase the tunnel splitting and ac- 
cess a regime in which the dynamics is expected to be 
coherent 0, Q ■ It is straightforward to show that the 
tunnel splitting depends on the transverse field to a high 
power, ~ D(Ht/ H a ) 2S f|]. This enables large variations 
in the magnitude of the tunnel splitting for small changes 
in the transverse field. Fig. ID shows the z-component 
of magnetization for symmetric and antisymmetric states 
as a function of the longitudinal magnetic field. Photon 
induced transitions (PIT) between symmetric and anti- 
symmetric states produce changes of the z component 
of magnetization when the longitudinal field is non-zero. 
When the microwave field, H ac , is parallel to the easy 
axis the transition rate is given by the following expres- 
sion: 

r = |(^ff oc ) 2 |(5|5 2 |A>| 2 /M (4) 

where (5*152 174) is the matrix element coupling the sym- 
metric and antisymmetric states, /(w) is a Lorentzian 
which characterizes the linewidth of the resonance: 

1 T 2 

7T 1 + (fj — LUO) T£ 

where u>o — E/h, and E is the energy level separation. 
T2 is the transverse relaxation time or decoherence time, 
which sets the width of the resonance. Note that the 
transition matrix element is maximum at zero magnetic 
field, when the states are symmetric and antisymmetric 
superpositions of \up) and \down) states. The matrix 



element decreases with longitudinal magnetic field and 
approaches zero when 2g^sH >> A, and states are sim- 
ple up and down spin projections. 

We have conducted experiments in a low temperature 
limit, ksT < E, in which there is a significant differ- 
ence in the population between the two lowest lying spin 
levels. A high sensitivity micro-Hall magnetometer is 
used to measure the z component of magnetization of a 
millimeter-sized single crystal of Ni 4 that is placed with 
one of its faces parallel to the plane of the Hall-sensor 
[Iril | (see fig. IE). The experiments were conducted at 
0.38 K in a He 3 cryostat that incorporates a 3D vector 
superconducting magnet, in which magnetic fields can be 
applied at arbitrary directions with respect to the axes 
of a crystal. A thin circular superconducting loop (</> = 2 
mm) is placed with its plane perpendicular to the easy 
anisotropy axis of the crystal (see Fig. IE). This loop 
shorts the end of a 2.4 mm coaxial line. A 50 GHz Vec- 
tor Network Analyzer is used to apply microwave radi- 
ation and enables characterization of the resonances of 
the loop. All the measurements were performed at loop 
resonances to maximize the power that is converted into 
ac magnetic field at the sample position. Further, pulsed 
microwave experiments were conducted using a pulse pat- 
tern generator to gate the microwave source. 

Figure 2A shows a magnetization curve recorded in the 
presence of continuous- wave (CW) radiation at 39.8 GHz 
while a transverse field of 3.2 T was applied. Peaks and 
dips are observed at opposite polarities of the longitudi- 
nal field demonstrating PITs between magnetic states of 
the molecules with opposite spin-projections. Note that 
at PITs the sample temperature increases by less than 
0.01 K, showing that these features are not due to sample 
heating. A careful inspection of the peaks shows a more 
complex structure: each peak is formed by two peaks, la- 
beled A and B. Measurements carried out at lower sweep 
rates show that peak B is also composed of two peaks, 
Bi and B2 (see Figs. 3 and 4). This is in agreement with 
longitudinal field EPR experiments, where different ab- 
sorption lines have been found and ascribed to molecules 
with slig htly different values of the axial anisotropy pa- 
rameter The inset of Fig. 2 A shows the difference 
between the magnetization measured with and without 
microwave radiation, M — M eq , versus field, for frequen- 
cies from 18 to 50 GHz and H T = 3.2 T. 

Classically the energy of these states would depend 
linearly on magnetic field and simply cross at zero field. 
However, the position of the peaks as a function of fre- 
quency does not depend linearly on magnetic field, par- 
ticularly near zero longitudinal field. This is evident in 
Fig. 2B, which shows the position of peak A versus lon- 
gitudinal field at different frequencies and as a function 
of the transverse field (2.2 to 3.6 T). The curvature of the 
energy splitting as the longitudinal field goes to zero is ev- 
idence of level repulsion and, with the measured magneti- 
zation changes, is a clear signature of quantum superposi- 
tion states with opposite magnetizations. The solid lines 
are fits of the data by direct diagonalization of the Hamil- 
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FIG. 2: (Color on-line) A) Magnetization versus longitudi- 
nal field in the presence of a transverse field, Ht = 3.2 T, 
while CW radiation of 39.8 GHz was applied. The inset shows 
M — M eq versus longitudinal field at different microwave fre- 
quencies. B) Field positions of the PITs of species A for 
several microwave frequencies and transverse fields. 

tonian of Eq. (1) using D A = 0.765 K, B = 7.9 x 10" 3 
K, C = 3.25 x lCT 5 K, g z = 2.3, g x = g v = 2.23 [H- 
The behavior of the B-peaks can be fit with parameters 
D B i = 0.735 K and D B2 = 0.745 K (not shown). These 
values are in excellent agreement with those determined 
from high frequency EPR experiments [TT| . 

In order to study the dynamics of the magnetization 
in the presence of radiation, we have conducted exper- 
iments with pulsed microwave radiation, in which the 
pulse time is of the order of the longitudinal relaxation 
time (~sec). We sweep the applied longitudinal mag- 
netic field at a rate of 8.3 x 10~ 5 T/s with H T = 3.2 
T around the position of a PIT. During this field sweep 
we apply pulses of 39.4 GHz with a 50% duty cycle and 
pulse times from 1 to 50 s. The results are shown in 
Fig. 3A. For the smallest sweep rate used in this ex- 
periment, 8 x 10~ 5 T/s, the field-width of the pulse is 
only 0.004 T (comparable to the size of the dots). This 
indicates that inhomogeneous broadening is not a signif- 
icant component of the lincwidth. The ratio between the 
magnetization when the radiation is OFF (solid circles) 
and ON (open circles) depends exponentially on the pulse 
width (Fig. 3B), M ff/M n = exp(-r/Ti). A fit to 
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FIG. 3: (Color on-line) A) Pulsed microwave results: 50% 
duty cycle pulses of 39.4 GHz radiation are applied while 
sweeping the longitudinal field across a PIT. The transverse 
field is 3.2 T. The magnetization was recorded while the pulse 
is ON (open circles) and OFF (close circles), for pulses widths 
from 1 to 50 seconds. The black curves were measured with 
99.99% (open black stars) and 0.01% (close black stars) duty 
cycles. B) The ratio between Mo ff and Mon for PIT A ver- 
sus pulse width for four different frequencies. The lines are 
fits to an exponential decay with a characteristic time n . C) 
n as a function of the longitudinal field. A decrease of the 
relaxation time is observed near zero field. 

this expression gives Ti, the longitudinal (or energy) re- 
laxation time, which governs the spontaneous decay of 
the excited state population. We have repeated this ex- 
periment for frequencies between 24.7 and 43.2 GHz at 
the same transverse field. The corresponding relaxation 
curves are shown in Fig. 3B. t\ is plotted as a function 
of longitudinal field in Fig. 3C. Remarkably, the relax- 
ation time is long and increases from 8 to 20 s as the field 
and thus frequency increase. This is contrary to general 
ideas that the relaxation time should decrease with fre- 
quency, because of the increasing phase space available 
for phonon generation [l3L Il4j . 

We have measured the change in magnetic response 
as a function of the power of the incident CW radiation 
(39.4 GHz) to estimate the decoherence rate. A trans- 
verse magnetic field of 3.2 T is applied and a longitu- 
dinal field is swept at a rate of 1.6 x 10 -4 T/s. The 
microwave power at the source was varied from -20 dBm 
to 2.5 dBm. The normalized change in magnetization 
[M — M eq )/M eq is shown as a function of the longitudinal 
field in Fig. 4A for different microwave powers. Fig 4B 
shows the amplitude for peak A versus power. The maxi- 
mum power applied to the sample produces transitions in 
20% of the molecules [3]. Larger magnetization changes 
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FIG. 4: (Color on-line) A) (M-M eq )/M eq for different powers 
with CW 39.4 GHz radiation applied to the sample. The 
longitudinal field was swept at f .6 x 10 -4 T/s and Ht = 3.2 
T. Note that the curves are offset for clarity. The red lines 
are fits to a Lorentzian with T2 = 0.6 ns. Figures B and C 
show the dependence of the amplitude and half width of PIT 
A, respectively. 

(up to 30%) have been obtained with higher transverse 
fields and at lower microwave frequencies (not shown) . If 
we consider that this amplitude corresponds to one of the 
peaks, which represents only a fraction of the molecules 
in the crystal, then the amplitude for this species is con- 
siderably larger. As an estimate, if each of the species 
represented one third of the total, then a net amplitude 
of 20% corresponds to 60% of molecules in that species. 
The amplitude versus power for peak A is shown in Fig. 
4B (assuming a 33% population). For this curve we can 
estimate the power necessary to saturate the resonance. 
We find 50 mW (source power) in this case (and about 
5 mW with H T = 3.4 T and / - 30 GHz, not shown). 
We have fit the results shown in Fig. 4A to a Lorentzian 



function f{ui) for the PIT A, as this feature is easily dis- 
tinguished from the others. We find that = 0.6 ns for 
the curve measured with the highest power (2.5 dBm) or 
a resonance quality factor of ~ 10. 

In Fig. 4C we show the dependence of PIT As half 
width at half maximum (HWHW) versus the microwave 
power at the source. The minimum width (lowest power) 
allows us to estimate a lower bound for the transverse re- 
laxation time. The results from the fit of the -12 dBm 
curve are shown in Fig. 4 A and give T2 > 0.8 ns. The ho- 
mogeneous linewidth represented by t-i will be less than 
that associated with inhomogeneous broadening. How- 
ever, we note that we recover the full amplitude of the 
peaks and dips during the pulsed-radiation experiments 
(see Fig. 3A), for sufficiently long pulses (r > 20 s). This 
indicates that ti ~ r^. A very similar decohcrence time, 
T0 = 1 ns, has been recently reported in high frequency 
EPR experiments carried out in a single crystal of a Mn4 
SMM dimer (S = 9/2) 

This experiment demonstrates the vast difference in 
time scales between the longitudinal and transverse relax- 
ations times of high spin superposition states. Clearly, for 
applications in quantum computing decoherence is a ma- 
jor challenge, which is likely amenable to synthetic strate- 
gies of nuclear isotope purification and further reduction 
in intermolecular interactions. Antiferromagnet molecu- 
lar magnets with an uncompensated electronic or nuclear 
spins that tracks the Neel vector for readout are also ex- 
cellent candidates for long coherence times in mesoscopic 
spin-systems and for use in quantum computation [17] . 
These will be the subject of future investigations along 
with studies of Rabi oscillations in molecular nanomag- 
nets induced by shorter and higher amplitude pulsed mi- 
crowave magnetic fields. 
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